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ABSTRACT 

We present GALEX near-UV {NUV) and 2MASS J band photometry for red se- 
quence galaxies in local clusters. We define quiescent samples according to a strict 
emission threshold, removing galaxies with very recent star formation. We analyse the 
NUV-J colour-magnitude relation (CMR) and find that the intrinsic scatter is an 
order of magnitude larger than for the analogous optical CMR (~0.35 rather than 
0.05 mag), in agreement with previous studies. Comparing the NUV-J colours with 
spectroscopically-derived stellar population parameters, we find a strong (> 5.5(t) cor- 
relation with metallicity, only a marginal trend with age, and no correlation with the 
a/Fe ratio. We explore the origin of the large scatter and conclude that neither aper- 
ture effects nor the UV upturn phenomenon contribute significantly. We show that the 
scatter could be attributed to simple 'frosting' by either a young or a low metallicity 
subpopulation. 
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1 INTRODUCTION 

The optical colour-magnitude relation (CMR) shows 
that brighter early-ty p e ga laxies are also redder 
ijVisvanathan fc Sandagg 1 19771 ). and is traditionally 
regarded as arising from the mass-metal licity sequence (of. 
lDresslej[l983 . iKodama fc Arimotd[l997l ). Gas loss, caused 
by supernova wind, occurs earlier in less massive galaxies. 
Therefore, a smaller fraction of gas is pro cessed before being 
expel led from the less massive galaxy (|Mathews fc Bakei 
1971 ) , resulting in lower average metallicities (|Larson 



19741 ). iBower. Lucev fc EUisI (|l992l ) found a very smaU 
intrinsic scatter in the U-V CMR (~0.05 mag) and, due 
to the sensitivity of the U band to the presence of young 
stars, interpreted this as a small age dispersion. Age 
and metallicity are observed to have the same effect on 
broadband optical colours, whereas spec tral line indice s 
can be used to break the degeneracy l|Worthevl 1 1994 ). 
iKuntschner fc Daviej (figgi) claimed that the CMR is 
driven by metallic ity variations with luminosity, although 
iNelan et all l|2005l ) found evidence for a strong age-mass 
relation in addition to this metallicity-mass trend (see also 
ICaldwell. Rose fc Concannonll2003l . iThomas et al.ll2005h . 
The ultraviolet-optical CMR for non-star-forming 



E-mail: t.d.rawle@dur.ac.uk 



galaxies has an intrinsic scatter an order of magnitude larger 
than its o ptical cou nterpart; ~0.5 mag compared to 0.05 
mag (e.g. lYi et all 2005). Hot young stellar populations 
dominate the ultraviolet (UV) flux for ~100 Myr after an 
episode of star formation (ten times longer tha n Ha emis- 
sion after star formation; iLeitherer et al. [l99^. The large 
intrinsic scatter in the UV CMR is therefore often inter- 
preted as differi ng quantities of very recent, albeit low level, 
star formation (|Ferreras fc Silk|[2000l '). 

In intermediate age populations (~l-3 Gyr), the near- 
UV (NUV; 2000-3000 A) flux is dominated by h ot stars on 
the main sequence turn-off (e.g. IO'Connel]|[l999l ). The sen- 
sitivity of the turn-off to the epoch of formation empha- 
sises the importance of the UV ban ds for age determination 
iDorman. O'ConneU fc Roodll2003l ). 

Old (~10 Gyr) metal-poor populations have a signifi- 
cant UV flux contribution from very hot (TeS ~ lOOOOK) 
blue horiz ontal branch (BHB) stars (|Maraston fc Thomas! 
2000; Lcc. Lee fc Gibsoru ^2002). However these tend to reside 
in globular clusters or galactic haloes (where Fe/H < — 1 ), 
where they are useful age indicators (|Kavirai et al.ll2007l ). 
rather than in relatively metal-rich elliptical galaxies. 

The UV picture is further complicated by the pres- 
ence of the ultraviolet upturn (or UV exc ess, UVX) phe- 
nomenon. First observed by ICodd (|l969l ). this unantic- 
ipated upturn dominates the far-UV (FUV; < 2000A) 
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in UVX galaxies. In contrast, the NUV can be decom- 
posed into two separate components: the blue end of the 
main sequence /subgia nt branch , and the UVX contribution 
ifDormanlQQTl .Bur stein et all (fl988) found that the UVX 
can sometimes be appreciable at wavelengths as long as 
2700A: for example, in NGC 4649 ~75 per cent of the NUV 
flux can be attributed to the UVX component. However, 
the UVX cannot be explained by the BHB population, as 
the temperature required to fit the upturn would be T^s > 
20000K, whereas BHBs are usually no hotter than Tcff ~ 
12000K (|0'Connelllll999l) . 

iBurstein et al.l further reported that FUV flux (as- 
sumed to trace the UVX) is strongly correlated with the 
Mg2 line strength (~metallicity) and also with the velocity 
dispersion, which i s a proxy for gala xy mass. However, more 
recent studies fe.g. lRich et al.ll2005^ have weakened the case 
for a strong UVX vs me tallicity relation. From analysis of 
internal colour gradients, lO'ConneU eiT^ |l992l l concluded 
that the FUV flux in most early-types originates from old 
stellar components. Drawing on these results, the source of 
the UVX is tentatively identified as hot, low mass, helium 
burning stars, such as extreme horizontal branch (EHB) 
or ' failed' AGB (AGB-manque) stars and the ir progeny 
(see l Yi. Demarque fc Oemlerlll997l . or the review lO'Connelll 
ll999F r 

T he Galaxy Evolution Explorer (GALEX; l aunched in 
2003; iMartin et al] |2005| . iMorrissev et all 120071 ) is revolu- 
tionising UV astronomy, with high resolution imaging in 
two bands: near- uhraviolet {NUV; \eff = 2310 A) and far- 
ultraviolet {FUV; Xeff — 1530 A). Usi ng analysis of both 
NUV-V and FUV-V vs B-V relations, lOonas et all (|2006l ) 
suggest that the FUV-NUV colour reflects an extension of 
the colour-metallicity relation into the UV, as well as de- 
ducing that ~10 per cent of elliptic als have residual star 
formation. Using the NUV- r colour, iKavirai et al.l (|2006l ) 
also flnd non-negligible young stellar p opulations in mor - 
phologically selected early-type galaxies. ISalim et al.l (|2007l ) 
derive star formation rates (SFRs) from broadband photom- 
etry dominated by the UV, and flnd good agreement with 
SFRs deduced from spectroscopic indices (predominantly 
using Ha). However, they also confirm that some galaxies 
with no Ha emission show signs of star formation in the UV 
bands and attribute this to post-starburst galaxies. 

Here, we build upon these previous studies by exploring 
the relationship between the NUV-J colour and spectro- 
scopic stellar population indicators for a sample of quiescent, 
red sequence galaxies in nearby clusters. This paper is organ- 
ised as follows. Section[2]describes our two red sequence sam- 
ples and associated 2MASS and GALEX datasets. The cri- 
teria used to remove galaxies with emission are described. In 
Section[3]we show that a large intrinsic scatter is found in the 
NUV-J colours of these quiescent cluster galaxies. Metal- 
licity is shown to be strongly correlated with the NUV-J 
colour, although there remains a large residual scatter. Sec- 
tion [4] discusses possible explanations for this scatter, show- 
ing that morphological abnormalities, aperture bias and the 
UV upturn do not contribute significantly. We investigate 
simple 'frosting' models with a low mass fraction of younger 
stars (or alternatively a low mass fraction population of low 
metallicity, blue horizontal branch stars), and show that 
these could account for the scatter. The uncertainties in the 



NUV K-correction are also discussed. Our conclusions are 
given in Section [5] 



2 DATA 

We use two complementary samples of red sequence galax- 
ies in local clusters: the first is explicitly red sequence se- 
lected by optical colour, and is a large sample, containing 
~10 times the number of galaxies; the second has the advan- 
tage of higher quality spectroscopy, and uses an emission line 
cut (Sec. 12. 3p to ensure a red sequence sample. ISmith et al.l 
12007) Fig. 1 demonstrates that Ha selection efficiently re- 
moves all galaxies bluer than the red sequence and is more 
restrictive than a cut on colour. 

2.1 Galaxy samples 

2.1.1 NFPS sample 

The NOAO Fundamenta l Plane Survey (NFPS: [SmLth et all 
|2004 iNelan et all |200i) is a study of X-ray selected clus- 
ters distributed over the whole sky and at redshifts between 
0.015 < z < 0.072. More than 4500 galaxies lying within 1 
Mpc of the centre of each cluster, and withi n 0.2 mag of 
the c luster red sequence on the B-R CMR (see ISmith et al.l 
|2004| Fig. 3), were observed spectroscopically. Of these, 3514 
have redshift, velocity dispersion and spectral line strength 
measurements (from 2 arcsec diameter fibres). 

2.1.2 Shapley Supercluster (SSG) sample 

The second sample of galaxies concentrates on the core of 
the Shapley supercluster (SSC; Abell clusters A3556, A3558, 
A3562 at z « 0.049). This sample consists of 541 galaxies 
selected from N FPS imaging but to a deeper limit {R < 18; 
ISmith et al.ll2007., l . Follow-up spectroscopy for these targets 
were obtained using 2 arcsec diameter fibres, equating to 2 
kpc at the distance of Shapley. A set o f three non-redundant 
line indices were fit to the models of iThomas et al.] (|2003l . 
|2004) in order to estimate age, metallicity (Z/H) and a- 
abundance (a/Fe) for each galaxy. The primary tracer of 
age in this scheme is H7F; for metallicity Fe5015 is used; 
Mgb5177 is the a-abu ndance indicator. T his method is de- 
scribed in detail iu iSmith et al.l l|in prep.h . 

2.2 GALEX and 2MASS data 

Galaxy Evolution Explorer {GALEX) near-ultraviolet 
{NUV) band images are available for 26 (from a total of 93) 
NFPS clusters. Due to a detector fault, only some clusters 
have associated far-ultraviolet {FUV) band images. Most 
of the images are from our guest investigator snapshot pro- 
gramme ( GALEX GI1_004) which targeted a subset of NFPS 
clusters with low ga lactic extinct i on an d having a large num- 
ber of spectra from lSmith et all l|2004h . In addition, a small 
number of images of comparable depth from the GR2 GI 
archive, medium imaging survey (MIS) and near galaxy sur- 
vey (NGS) have been used. Table[T]lists all the images anal- 
ysed along with their centre position and stacked exposure 
time. GALEX images have a 1.25° diameter, but only the 
central 1.2° field has been analysed due to the poor image 
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Table 1. GALEX NUV images, cz (kms ^) is cluster redshift in the CMB frame. A single exposure is used unless column six indicates 
the number of coadded images. Total exposure time for coadded images is given for NUV and (where the image is available) FUV. 



GALEX image Centre RA Centre Dec Cluster(s) czqmb Images tcxp (sees) 

(J2000) (J2000) (kms-i) coadded NUV FUV 



GI1_004001_A2734 


00 


11 


21 


6 


-28 


51 


00 


A2734 


18249 


3 


3535 




GI1_067001.UGC0568.0003 


00 


55 


08 


9 


-01 


02 


47 


A0119 


12958 




1556 


3024 


MISDRl.16976.0422 


01 


14 


29 


3 


+15 


01 


46 


A0160A 


12794 




1444 


1444 


GI1.004002.A0262.0001 


01 


52 


45 


6 


+36 


08 


58 


A0262 


4464 




1698 


1698 


NGA.NGC1058 


02 


43 


26 


6 


+36 


25 


39 


A0376 


14371 




1265 


1265 


GI1.004003.A3104.0001 


03 


14 


21 


6 


-45 


25 


12 


A3104 


21560 




1588 


1588 


GI1.004027.A3158.0001 


03 


42 


57 


6 


-53 


37 


48 


A3158 


17542 




1026 


887 


GI1.004004.A3266 


04 


31 


24 





-61 


26 


24 


A3266 


17713 


2 


1441 


1441 


GI1.004005.A0548 


05 


46 


40 





-25 


37 


21 


A0548A/B 


12439 


4 


3155 




GI1.004006.A3376 


06 


01 


43 


2 


-39 


59 


24 


A3376 


14016 


2 


2176 


2176 


GI1.004007.A3389.0002 


06 


21 


57 


7 


-64 


57 


35 


A3389 


8075 




976 




MISDRl.24335.0270 


10 


13 


07 


7 


-00 


23 


25 


A0957 


13849 




1703 


1703 


GI1.004025.A3528.0001 


12 


53 


57 


6 


-29 


13 


48 


A3528A/B 


16764 




1697 




GI1.004008.A1644.0002 


12 


57 


12 





-17 


24 


36 


A1644 


14478 




1163 




GIl.009003.HPJ1321m31.0001 


13 


21 


05 


8 


-31 


32 


20 


A3556 


14660 




1615 




GI1.004010.A3556 


13 


25 


26 


1 


-31 


36 


07 


A3556, A3558 


14660 


2 


1805 




GI1.004011.A3558.0001 


13 


27 


57 


6 


-31 


30 


00 


A3558 


14660 




1676 




MISDRl.33707.0586 


14 


42 


46 


3 


+03 


39 


11 


MKW8 


8449 




1698 


1698 


GI1.004016.A1991.0001 


14 


54 


31 


1 


+18 


38 


23 


A1991 


17741 




967 




GI1.004026.A2063.0001 


15 


23 


36 





+08 


36 


34 


A2063 


10444 




1512 




NGA.NGC6166 


16 


28 


39 


9 


+39 


33 


24 


A2199 


8872 




1437 


1437 


GI1.004020.A3716 


20 


51 


57 


5 


-52 


46 


48 


A3716 


13141 


2 


3179 


1689 


GI1.004021.A2399.0004 


21 


57 


19 


1 


-07 


47 


59 


A2399 


17046 




1322 




GI1.004022.A2589 


23 


23 


57 


6 


+16 


46 


47 


A2589 


12001 


3 


4252 




GI1.004023.A4059.0003 


23 


56 


59 


9 


-34 


45 


35 


A4059 


14660 




635 





quality at the edges. GALEX images have a plate scale of 
1.5 arcsec pixel" 1 and a PSF FWHM of ~5 arcsec. 

We analyse infrared tiles from the J band (Ae// = 
1.25^m) of the Tw o Micron AU Sky Survey (2MASS; 
ISkrutskie et al.ll2006l ). We measure directly from the tiles, 
rather than adopt pho tometry from the 2 MASS extended 
source catalogue (XSC; Ijarrett et al1l200d ). as some of our 
target objects are unresolved. Additionally, we require the 
PSF to match that of the NUV images (2MASS J PSF 
FWHM ~3 arcsec) and therefore Gaussian smooth the 
2MASS tiles before analysis. 2MASS tiles have a plate scale 
of 1 arcsec pixel" ^. 

SExTRACTOR (|Bertin fc Arnout jliggd ) was employed 
in dual image mode to detect all objects in the J band. Pho- 
tometry was measured for all sample targets in both bands 
using a range of matched apertures: a Kron-type aperture 
(SExtractor's MAG.AUTO) and seven apertures 3-34.5 
arcsec in diameter (MAGjVPER). Throughout this work, 
J band Kron apertures are used for total luminosity, and 
matched 12 arcsec diameter apertures for NUV-J colours. 
Our J band photometry is in good agreement with 2MASS 
(~0.2 mag RMS) for objects in the XSC. 

Targets with SExtractor apertures flagged as trun- 
cated, or with a deblending error, have been removed from 
the sample, and only confirmed cluster members with red - 
sh ift and log a measu rements within the lSmith et al.l (|2004l ) 
or lSmith et all (|2007l ') datasets (NFPS and SSC respectively) 
are used in the analysis. Table[2]lists the number of galaxies 
in the samples at this stage, and after subsequent restric- 
tions. 

All colours and magnitudes are measured in the 



Table 2. Size of the galaxy samples, following the selection cri- 
teria applied. 





NFPS 


SSC 


Original sample 


4527 


541 


...with usable NUV and J photometry 


1493 


307 


...with cz and logcr data 


990 


267 


...after emission line cut 


920 


156 


after optical apparent magnitude cut^ 


544 


101 


with stellar population parameters^ 




87 


with FUV and H7F data^ 


222 





NOTE - see ^(SeclSlTJ, ^(SeclsSll or ^(Sec li^jl 



AB system (Jab = Jvcga + 0.91; iBlanton et al] 120051 ) . 
and have been corrected for galactic extinction using 
the reddening maps of ISchleeel. Finkbeiner fc DaviesI 
(1998); Afuv=8.29 xE(B-V), Anuv=8.87xE ( B-V) 
Aj=0.902x E{B-V) (jCardelh. Clavton fc Math5 Il989l : 
Schlegel et al.). 

Estimates of the K-correction in the UV bands are cur- 
rently derived empirically from poorly constr ained spectra, 
and t herefore subject to large uncertainties. iKavirai et al.l 
(|2006l . fig. 22) estimate that the NUV-i correction would 
be ~0.1 mag throughout the range in redshift considered 
here. We do not apply a K-correction in this study. This 
issue is addressed further in Section [4.41 
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Figure 1. Colour-magnitude diagram for tiie samples (NFPS 
upper panel; SSC lower panel). Grey points are those removed by 
emission line criteria detailed in Sec. 12.31 The larger tail of grey 
points for SSC results from the lack of an explicit colour cut on 
the original sample. Solid lines are GALEX detection limits and 
dashed lines are the applied J band apparent magnitude cuts. 
Median error bars are ~0.05 mag in J band and ^^0.1 mag in the 
colour. 



2.3 Emission-line cuts 

In order to construct a sample of emission-free red sequence 
galaxies, a restriction is made on emission line strengths. 
Mindful of the effect of nebular emission 'fill in' for the age- 
sensitive Balmer lines (H/3, H7 and H(5), the preferred cut 
is on the Ha line. Unfortunately, Ha was not measured for 
most NFPS galaxies, so the selection criteria of the origi- 



nal NFPS reduction has been adopted (|Nelan et al.l 12005 
Specifically, this involves a cut on the H/3 equivalent width, 
EWiYi-P) > —0.6 A (negative i?W denotes spectral line emis- 
sion) supplemented by a cut on OHI A5007, SW(OHI A5007) 
> —0.8 A. The cut on the SSC sample, which has Ha mea - 
surements, follows the prescription of 'S mith et al.l (|2007l ). 
which uses ^"^(Ha) > —0.5 A (approximately equivalent to 
EWiUP) > -0.2 A). The OHI cut would not remove any 
additional galaxies. These cuts ensure red sequence subsam- 
ples, free of galaxies with a sizeable star formation compo- 
nent or an optically strong active galactic nucleus. The data 
and photometry for the resulting subsamples are reported 
in Tables E] and H 



3 RESULTS 

3.1 NUV—J colour relations 

Figure [1] shows the NUV-J colour-(apparent) magnitude 
diagram for the NFPS and SSC samples. Shown in grey are 



the galaxies removed by the emission cut from Section [2.31 
Prior to the emission cut, SSC has a larger colour range than 
NFPS. This is because NFPS was explicitly selected on B-R 
colour while SSC only on total R band apparent magnitude. 
The Ha cut efficiently removes the very blue objects. 

All target galaxies were within the 2MASS J band de- 
tection limit, but the NUV band has a 5a detection limit of 
22.5-23.5 mag, depending on the coadded image exposure 
time. Figure[T] (upper panel) shows that for the NFPS, these 
limits result in a bias against faint red galaxies. Assuming 
the brightest NUV detection limit, a sample cut is applied 
to the J band apparent magnitude at 15.0 mag. For the SSC 
sample, all of the targets appear on just two NUV images 
with similar exposure times, and a {5a) detection limit of 
23.2 mag. Although in practice only two SSC targets have 
a non-positive fiux in the NUV band (compared to ~8 per 
cent in NFPS), for consistency, SSC has been treated in a 
similar manner, with a cut applied at 15.8 mag. These cuts 
are shown as dashed lines in Figure [1] and ensure a complete 
sampling of the colour range over the selected luminosity 
interval. 

The colour-(absolute) magnitude diagrams for the low- 
emission galaxies are presented in Figure (2] There is a corre- 
lation between the luminosity and colour; brighter galaxies 
tend to be redder. However, there is a large scatter; rms dis- 
persions of 0.37 and 0.30 mag for NFPS and SSC samples 
respectively. The smaller scatter within the SSC sample is 
probably due to the slightly more restrictive Balmer emis- 
sion line criteria (see Sec. 12. 3p . The scatter in each sample 
does not increase by more than ~10 per cent unless the 
cut criteria are relaxed beyond an equivalent width of -1 
A. Only 5 per cent of the scatter can be accounted for by 
photometric measurement error. As intrinsic scatter domi- 
nates, all correlations in this study are computed without 
error weighting. 

Table [S] summarises the CMRs for our two sam- 
ples. Additionally, the C MRs measured in three 
previous studies (|Yi et aL 



iHaines. Gargiulo &: Merluzzil 



20051 . iBoseUi et al.] |2005| . 



20071 ) are shown for compari- 



son. Di fferent sam ple selections were used for each of these 
studies. lYi et al] use a 'UV-weak' early-type galaxy sample, 
selected from SDSS by concentration index and luminosity 
profiles, and then by the fiux ratios F{NUV)/F{r) and 
F{FUV)/F(r) both being less than 0.07. The sample 
covers a J b and luminosity range comparable to our work. 
iBoselli et alj use a volume-limited sample of galaxies in 
the Virgo cluster, with a subsample defined as elliptical 
by visual classification. Haines et al. use a volume limited 
sample of local galaxies from SDSS, with the subsample 
(labelled 'passive red sequence galaxies') restricted by the 
emission line criteria EW{}la) > —2 A. 

Despite these selection differences, our derived relations 
are in good agreement with the previous studies. The scat- 
ter is consistent (~0.3-0.5 mag) given the different sample 
definitions, and considerably large in comparison to that of 
the optical OMR (~0.05 mag; Bower et al. 1992). 

Velocity dispersion provides an alternative mass proxy 
to luminosity, and, for opt ical colours, the a correlation 
appears more fundamental l|Bernardi et al J [20051 ). Figure [3] 
presents the samples in terms of their NUV-J colour and 
logo". There is a clear correlation, with slopes of 0.73 ± 0.11 
and 0.65 ±0.17 for NFPS and SSC respectively, but the rms 
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Table 3. Data for galaxies in the NFPS sample. Galaxy position is encoded in the ID. czcMB (kms~^) is the mean cluster redshift in 
the CMB frame, czi^^i (kms~^) is galaxy redshift in the heliocentric frame. Magnitudes (Kron-type apertures) and colours (12 arcsec 
diameter apertures) are in the AB system and have been galactic extinction corrected, but not K-corrected. 



Galaxy ID 




cluster 


galaxy 


apparent J 


NUV-J 


FUV-J 






C^CMB 


C2hcl 








NFPJ043305.0- 


-612235 


17714 


16144 


15.750 ± 0.071 


6.732 ± 0.142 


7.205 ± 0.250 


NFPJ043306.7- 


-612614 


17714 


17702 


14.898 ± 0.047 


7.091 ± 0.102 


7.497 ± 0.184 


NFPJ043307.6- 


-611338 


17714 


16356 


14.659 ± 0.042 


6.979 ± 0.095 


7.659 ± 0.195 


NFPJ054415.7- 


-255429 


12939 


10674 


14.693 ± 0.045 


6.886 ± 0.077 




NFPJ054431.6- 


-255550 


12939 


13247 


14.783 ± 0.047 


6.767 ± 0.086 





Full content of this table is available in the electronic version 



Table 4. Data for galaxies in the SSC sample. As in Table|3] Supcrcluster mean c^cMB 
parameters given where available (age in Gyr). 



14660 km s for all galaxies. Stellar population 



Galaxy ID 




apparent J 


NUV-J 


log(age/Gyr) 


[Z/H] 


[a/Fe] 


NFPJ132418. 2-314229 


13948 


14.197 ± 0.035 


7.088 ± 0.085 


0.94 ± 0.03 


0.27 ± 0.02 


0.24 ± 0.02 


NFPJ132423. 0-313631 


14642 


14.847 ± 0.048 


6.706 ± 0.103 


0.89 ± 0.05 


0.11 ± 0.03 


0.15 ± 0.02 


NFPJ132425. 9-314117 


13888 


14.380 ± 0.038 


6.604 ± 0.079 








NFPJ132426. 5-315153 


14922 


14.676 ± 0.043 


6.575 ± 0.081 


0.95 ± 0.04 


0.23 ± 0.03 


0.24 ± 0.02 


2MASXJ13250387-3132449 


14266 


14.792 ± 0.049 


6.626 ± 0.109 









Full content of this table is available in the electronic version 
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Table 5. NUV-J CMRs. 





Absolute J (mag) 

Figure 2. NUV-J colour— absolute magnitude diagram for both 
samples, with best fits in black. R elations from previ ous stud- 
ies in grey:lYi et al.l l2005h=dashed,lBoselli et al.l (l2005l) =dotted. 
iHaines. Gargiulo &: Merluzzil |20o3)=dash-dotted. Median error 
bars are ~0.05 mag in J band and ~0.1 mag in NUV-J. 







original 




transformed 


rms 




X 


y 




a 


b 
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vironments than our samples. Their blue galaxy fractions 
are plotted in Fig. 2] and show a markedly higher fraction 
for a given velocity dispersion. For example, our two sam- 
ples have a blue fraction of ~40 per cent only at t he low- 
est sigma (logcr < 1.8), while the ISchawinski et al.l sample 
reaches this blue fraction at logcr — 2.2. At face value, this 
result suggests a large difference between field and cluster 
galaxy populations. However , the differences i n sample se- 
lection have to be considered (jSchawinski et ahl use a sample 
selected on morphology), which is beyond the scope of this 
study. 



scatter (0.36 and 0.32 mags) is indistinguishable from that 
of the CMR. 

Figure |4] shows the fraction of 'blue' galaxies as a func- 
tion of log a for our two samples. Blue galaxies are d efined 
by NUV-J < 6.4 mag, as used in ISchawinski et al.l (|2006l ) 
who study SDSS galaxies probing to much lower density en- 



3.2 Stellar population parameters in the Shapley 
sample 

In order to investigate the physical origin of the large in- 
trinsic scatter found in the NUV-J colour, we examine the 
relationship between colour and the stellar population pa- 
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Figure 3. NUV-J vs logo- for NFPS (upper panel) and SSC 
(lower panel). Median error bars are ~0.02 and ~0.01 for logtr in 
the two samples respectively and ^O.l in the colour. 




log <j 

Figure 4. Fraction of galaxies with NUV-J < 6.4 as a function 
of logfT. NFPS: squares/s olid line; SSC: tr iangles /dashed line. 
Estimates from Figure 1 of lSchawinski et al.l (HoOQ) are given for 
comparison (circles/dotted line) and show a much steeper increase 
in the number of blue galaxies with decreasing a. 



rameters (age, metallicity, a-abundance) for the SSC sam- 
ple. 

In the following analysis, involving only the SSC sam- 
ple, the derived limit on J band apparent magnitude (Sec. 
13. 1|) is not applied, allowing a larger number of galaxies. 
(It should be noted that the subsequent conclusions are ro- 
bust against the use of the J band cut.) The emission line 
selection criteria are retained, and the sample is further re- 
stricted by the overlap of the photometric and spectroscopic 
datasets (see Table [2]). 

NUV-J vs log(age) is shown in Figure [5] (upper panel). 
Galaxies with nebular emission have been removed via the 
emission line criteria. Therefore, if the UV sources also con- 
tribute to the optical flux and have strong Balmer lines, the 
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Figure 5. NUV-J (upper) and B—R (lower; from 
|m ercurio et al.h colours vs log(age) for the SSC sample. 
Symbol size represents J band luminosity; larger=brighter. Note 
that the lower panel has a much smaller range in the vertical 
scale. 



NUV-J and age for the remaining objects would be corre- 
lated. We find only a marginal correlation between age and 
NUV-J, with a slope of 0.46 ± 0.25 and an intrinsic scatter 
(after accounting for the measurement error) of 0.33 mag. 
From the upper panel of Fig. [5] it is also apparent that more 
luminous, and by inference larger, red sequence galaxies (J 
band luminosity is shown in Figs.[5H7]by the symbol size) are 
not solely confined to the redder NUV-J colours, although 
they do tend to be the oldest. 

Figure \S\ (lower panel) shows the relation between age 
and B-R colour for the same sample of galaxies (note the 
much smaller range in the vertical scale). There is a strong 
correlation; a slope of 0.21 ±0.04 and rms dispersion of 0.05 
mag. This figure confirms that the NUV scatter is not due 
to contamination by optically blue galaxies. 

The analogous NUV-J correlation with the metallic- 
ity is given in Figure [S] There is a strong trend between 
[Z/H] and NUV-J, with a slope of 1.27 ± 0.23 and an 
rms dispersion of 0.32 mag (~90 per cent of which is in- 
trinsic scatter). The lack of galaxies to the top left is not 
the result of a selection effect. Of the two target galaxies 
undetected in the NUV band (see Sec. I3.1|l . only one heis 
a low metallicity. Assuming the non-detection is due to a 
redder-than-average colour, this would add a single galaxy 
to the upper left of the plot, but would not significantly 
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[Z/H] 

Figure 6. NUV-J vs metallicity [Z/H]. Symbol size represents 
J band luminosity; larger=brighter. 




0.1 0.2 0.3 0.4 0.5 

[«/Fe] 

Figure 7. NUV-J vs a-abundance [o/Fe]. Symbol size repre- 
sents J band luminosity; largcr=brighter. 



affect the fit. In general, tlie most luminous galaxies form 
a ridgeline at redder colours, althougii tiiere are significant 
brigfit outliers to this trend (the most obvious being the 
bright, blue, metal-rich galaxy on the right; NFPJ132729.7- 
312325). Lower metallicity g alaxies tend to be blue r and 
less luminous. This supports iRampazzo et all l|2007l ). who 
use simulations to predict a correlation of NUV-IK colour 
with metallicity, but little dependence on age in popula- 
tions greater than 2-3 Gyr after a star formation episode. 
However, the slope of our observed metallicity trend is 2-3 
times weaker than that derived from theoretical spectra by 
iDorman. O'Connell fc R ood (200^). 

Model evoluti onary tracks suggest that stellar evolutio n 
depends on a/Fe (|Salasnich et al.li2000l : fPotter et al.ll2007^ . 
However, we find no discernible relation between the NUV- 
J colour and the a-abundance in this sample of red sequence 
galaxies (rms scatter of 0.37 mag; Fig. [T}. 

Figure |8] shows the residuals from the NUV vs metal- 
licity relation (Fig. |S} against log(age). Although the resid- 
uals are more strongly correlated with age (~3(t) than the 
colours themselves are, the rms dispersion is only reduced to 




Figure 8. NUV— J residuals from the [Z/H] relation vs log(age). 
Symbol size represents J band luminosity; larger=brighter. 

0.3 mag (~0.25 mag intrinsic scatter). This is in contrast to 
the case of the B-R colour where the majority of the scat- 
ter can be attributed to age and metallicity l|Smith et al.l 
lin prep.h . 



4 DISCUSSION 

This section explores possible causes of the large intrin- 
sic scatter observed in the NUV -J colour. We investigate 
aperture bias, the UV upturn phenomenon and 'frosting' by 
young or metal poor subpopulations. Additionally, we com- 
ment on the uncertainty in the NUV band K-corrections. 

4.1 Aperture bias/morphology 

The spectroscopy used to estimate the stellar parameters 
is derived from 2 arcsec diameter fibres, whilst the NUV- 
J colours are from 12 arcsec diameter aperture photometry. 
Eliminating aperture bias completely would require matched 
apertures for the photometry and spectroscopy, but unfor- 
tunately the PSF of the 2MASS and GALEX images is too 
large for reliable 2 arcsec aperture photometry. 

Broadband colours and spectroscopic measurements 
(e.g. iTamura fc Qhtal l2004l . ISanchez-Blizauez et al.l l2007h 
show that elliptical galaxies generally have flat radial pro- 
files in age, and regular metallicity gradients (decreasing 
[Z/H] with radius). The aperture effect therefore flattens the 
NUV- J vs Z/H relation, as larger, redder galaxies will tend 
to exhibit higher metallicities within the fibre. However, the 
effect is small, with only ~0.1 dex change in metallicity over 
a 1 dex difference in aperture radius. 

Inspection of the galaxies in high resolution images can 
ascertain whether there are morphological peculiarities, or 
neighbouring objects, contributing to an enhanced large ra- 
dius NUV flux. Figure |9] shows colour vs age and metallic- 
ity (as in Figs.[SHB}, highlighting the few galaxies that have 
high resolution HST- ACS images available. Fortunately, one 
of these objects is the most obvious outlier in the whole 
sample (NFPJ132729.7-312325); it is the reddest in B-R, 
but has an unusually blue NUV-J colour for such a metal- 
rich, luminous galaxy. In the HST- ACS image (Fig. llOp. the 
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log(age/Gyr) 

Figure 9. Colour vs stellar parameters, highlighting galaxies with 
available HST-ACS images. Symbol size reflects J band luminos- 
ity (larger=brighter) . In all cases, ACS confirms early- type mor- 
phology, uncontaminated by neighbours. 




Figure 10. HST-ACS image of the outlier NFPJ132729.7- 
312325, with a 12 arcsec diameter aperture marked. 



galaxy appears to be a large elliptical with no abnormalities. 
There are also no obvious contaminating objects that could 
be responsible for an anomalous blue colour. 

On close examination, the other targets with HST-ACS 
images also appear to be 'normal' ellipticals of various sizes, 
with no obvious peculiarities or large radius NUV contrib- 
utors. We conclude that a large scatter in NUV-J colour 
is present even in objects with confirmed regular early-type 
morphologies and no contaminants. Additionally, a scatter 
of ~0.3-0.35 mag is still obtained when smaller matched 
apertures (4.5 arcsec diameter) are used for the colour, de- 
spite the probability of a contaminating neighbour being 
reduced by ~85 per cent. 



4.2 UV upturn 

Another possible explanation for the scatter in NUV-J 
colo ur is NUV contamin ation by the UV upturn (UVX). 
The iThomas et all (|2003l ) models used in the stellar popu- 
lations parameter calculations include low-metallicity blue 
horizontal branch sub-populations and thermal pulsing 
asymptotic giant branch stars, but they do not include the 



low mass, metal-rich, helium burning stars with small en- 
velopes currently t hought to be th e most likely candidate 
for the UV upturn (|0'Connelllll999l '). The UV upturn is one 
of the most heterogeneous photometric properties of old stel- 
lar populations in early-type galaxies, with a s pread of up to 
~4 magnitudes in the FUV ( 0'Conneiilll999l ). so certainly 
seems a plausible explanation for the scatter. 

The hot UVX component appears in the spectra as a 
smooth continuum with an absence of emission and absorp- 
tion lines. Hence, the FUV flux from a galaxy can be used 
to estimate the exten t of the contribution in the NUV band 
IIBurstein et al .Il988f ). Unfortunately, none of the SSC galax- 
ies have available FUV photometry, so in order to estimate 
the extent of the UVX contribution, the corrections are cal- 
culated for all objects in the NFPS sample that have the 
necessary FUV data. However, stellar population ages have 
not been derived for individual NFPS galaxies, so instead 
we compare the corrected NUV-J colours to the traditional 
'age-tracing' Balmer line H7F. This line does not trace age 
cleanly, being affected to a small degree by the metallicity. 

H7F vs NUV-J for the NFPS sample is shown in Fig- 
ure 111! Given the relative strengths of the colour vs. age 
and [Z/H] relations (Figures [5] and [B]), the correlation seen 
here is most likely a reflection of the metallicity, rather 
than age, dependence of H7F. For comparison, the observa- 
tions have been overlaid by age / meta llicity grids constructed 
from the models of lMarastonl (|2005l . M05; upper panel) and 
iBruzual fc Charlotl l|2003l . BC03: lower panel). Both grids 
lie redward of the observed data, most likely due to nei- 
ther model including EHB stars, although the BC03 grid 
provide s the better description of the data. Previous stud- 
ies (e.g. ISalim et al.l[2007 'l have noted that while the BC03 
models do not explicitly include EHB stars, the UV light 
from old stellar populations (primarily post-AGB stars) re- 
produce several of the correlations found in observational 
data, includ ing the relation be tween FUV-NUV and B-V 



colo urs (see Donas et al.| 20061. particula rly fig. 5) 



iDorman. Q'Connell fc RoodI (|2003l ) introduced correc- 
tions to their NUV-V colours by considering the relative 
contributions of the hot and cool (i.e. non-UVX) compo- 
nents at FUV, NUV and V. They assumed a negligible 
contribution of the cool component to the FUV flux, but 
allowed hot stars to contribute at V. Here, we use analo- 
gous corrections for the NUV-J colours, which are simpler 
because it is safe to assume no UVX contribution at J. The 
corrected colours are given by 



{NUV - J)corr = -2.5 X Zoff[10-"-^(^^^-'''=^> 
-a X io-«-4(^^^-^'=-l 



(1) 



where a ~ 0.3 is the ratio of NUV to FUV flux for the 
hot component (appropriate for a T^h — 24000 K star, see 
Dorman et al.). 

The resulting corrections are shown as vertical lines in 
Figure [12] (where the points indicate the value of the cor- 
rected colour). The corrections move the colours redward as 
expected, but are only of the order of ~0.2 magnitudes, and 
do not reduce the scatter in the NUV-J (0.29 mag before 
corrections, 0.30 mag after). Obviously there are no extreme 
UVX galaxies with large hot component contributions to 
the NUV fiux in this sample, and as such the UVX phe- 
nomenon is unlikely to responsible for the scatter in NUV- 
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Figure 11. NUV-J (without K-correction) against H7F Balmer 
index for NFPS galaxies. Upper panel: M05 grid (age={5,12,15 
Gyr}, [Z/H]={-0.33,0.00,+0.35}). Lower panel: BC03 grid 
(age={5,12,15 Gyr}, [Fe/H]={-0.33,+0.09,+0.55}). (5,-0.33) is 
the bluest (age,[Z/H]) grid point. Median error bars are shown 
top right. 



J for NFPS. We speculate that the UVX effects in the SSC 
sample are smaller because, on average, the galaxies have 
lower luminosities. 



4.3 Stellar population 'frosting' 

The simple models can be generalised by constructing com- 
posite stellar populations. A common invocation of this is 
residual star formation in the form of 'fr osting' galaxies wit h 
a small mass fraction of young stars (|Trager et al.ll200d ). 
which manifests itself in the observables as a bluer colour 
and a younger age than the base population. Frosting should 
not affect the spectroscopic metallicity or a-abundance sig- 
nificantly as these arc primarily driven by the larger, older 
population (jSerra fc Trageij|2007l ). 

Figure [13] shows the extent to which frosting can ac- 
count for the scatter, after the effect of metallicity has been 
removed using the correlation of Fig. [S] Simple stellar popu- 
lation (SSP) tracks from both the M05 and BC03 models are 
shown. The slope of the BC03 track fits the observed red en- 
velope well, whereas M05 predicts a much steeper variation 
with age. Allowing for the systematic effects of aperture bias 
and the UV upturn (~0.3 mag bluer in observed NUV-J), 




H7F 

Figure 12. NUV-J (without K-corrcction) against H7F Balmer 
index for galaxies with available FUV data. Points indicate the 
corrected colour, with the 'tail' showing the correction vector. 
BC03 grid as in Fig. 1111 (lower panel). Median error bars are 
shown top right. 



~30 per cent of the galaxies cannot be accounted for by the 
SSP model. 

Vectors of the frosting effect on NUV-J colour and 
spectroscopic age have been calculated using the BC03 mod- 
els. Frosting of a 15 Gyr base population by a 1.5 Gyr 
population with a mass fraction /i=0.03, and by a 0.7 Gyr 
population of /i=0.01 are shown. The M05 models result in 
marginally steeper frosting vectors, as would be expected 
given the steeper SSP track. The vectors show that frosting, 
even at a modest level of 1-2 per cent for a ~1 Gyr popu- 
lation, could account for a sizeable portion of the scatter in 
the NUV-J colour. A young population of this size would 
not be apparent in the B-R colours. 

The spectroscopic age is sensitive to frosting via the in- 
creased hot-star contribution to the Balmer lines, and for a 
given change to the spectroscopic age, the UV colours are af- 
fected more strongly by frosting than by lowering the age of a 
single-burst population. This supports the assertion that UV 
colou rs are partly dependent on low level recen t star forma- 
tion (|Ferreras fc Silkll2000l . iKavirai et allbooel . [ Salim et al.l 
l2007h . However, there are spectroscopically old galaxies with 
blue colours (~10 per cent of the sample) which cannot be 
accounted for by the frosting scenario describe d above. Ad- 
ditionally, Rose Call index results ( Smith ct al. i n prep. I) ap- 
pear not to support the presence of young stellar populations 
(< 1.5 Gyr) in the majority of red sequence galaxies. 

For most of this study, we have neglected the blue hor- 
izontal branch (BIfB) as cluster red sequence galaxies have 
[Z/H] ~0. However, it is possible that a low metallicity pop- 
ulation with a BHB morphology m ay be present in some 
galaxies (|Maraston fc Thomas! [2OO0I ) . Therefore, an alter- 
native 'frosting' scenario consists of a low mass fraction, 
old, low metallicity, BHB stellar population embedded in a 
[Z/H]=0 galaxy. Estimates from the M05 models show that 
'frosting' by a 4-5 per cent mass fraction population with 
[Z/H]=-1.35 could give an NUV-J colour ~1 mag bluer, 
while changing the B-R by only ~0.03 mag. The effect of the 
low-metallicity frosting on the derived spectroscopic age is 
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Figure 13. Simple stellar population tracks from the M05 and 
BC03 models shown in the context of the NUV-J colour (cor- 
recting to [Z/H]=0 using the correlation from Fig.|6]l vs log(age). 
Vectors indicate two BC03 frosting scenarios with minor popu- 
lations of {M,age} = {0.03,1.5Gyr} and {0.01,0.7Gyr} (solid and 
dashed vector respectively). 



likely similar to the effect of frosting by a young component, 
since both cases are driven by increased A-star contribution 
to the Balmer lines. 



4.4 NUV K-corrections 

K-corrections for the UV bands are poorly constrained, as 
the spectral shape of galaxies at wavelengths shorter than 
^3000A is not well known. Previous studies of the UV CMR 
largely 

(120051, who estimate corrections of 0.1 



;ely ignore the K-correction, or in the case of lYi et al.l 
35l . who estimate corrections of 0.1 — 0.2 mag for z = 
— 0.25) apply correcions based on the luminosity distance 
wit hout considering the spectral shape of galaxies. 

iKavirai et aL I ||2006D compute NUV K -corrections for 
best fit model SEDs derived from SDSS and GALEX pho- 
tometry. Corrections of ~0.1 mag were found for redshifts z 

< 0.1. However, corrections for their model of a 9 Gyr old 
simple stellar population are much larger (0.4 — 1.0 mag for 
0.04 < z < 0.11). 

As a n illustration, we calculate K-corrections from the 
models of iMarastoiil (|2005l . M05) at various redshifts (0.015 

< z < 0.072), both for simple populations and galaxies in the 
'frosting' scenario described in Sec l4.3l K-corrections of 0.2- 
1.2 mag are obtained, depending on the metallicity and age 
(increasing either increases the correction), and on the mass 
fraction of the young stellar component. T he corrections are 
dominated by the 2640A spectral break (cf lEisenstein et"al] 
I2OO3I ). which i s redshifted comp letely out of the NUV band 
by 2 « 0.07. In lEisenstein et all , who study the average spec- 
tra of 726 luminous, red, SDSS galaxies at 0.47 < z < 0.55, 
the break appears less prominent than in the M05 spectra 
(a ~50 per cent drop in flux as opposed to ~80 per cent). 

Fig [14] shows the regression line residuals from the 
colour-magnitude diagram for the NFPS sample (Fig [2] up- 
per panel) plotted against the galaxy redshift in the he- 
liocentric frame. The solid line shows the expected trend 
if the M05 K-corrections were necessary, but not applied. 
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Figure 14. CMR residuals for the NFPS sample (Fig (2] upper 
panel) vs heliocentric galaxy redshift. The filled squares and solid 
line show the predicted trend if the K-corrections from the M05 
models were necessary, but not applied. 



This implies the uncorrected residuals should have a steep 
correlation with redshift, which is not obser ved. A flatter 
K-cor rection is preferr ed by the data, as in iKavirai et al.l 
(|2006| ). Along with the' Eisenstein et al.l spectrum described 
above, this result highlights the uncertainty in the UV K- 
corrections (as well as the problems of stellar population 
models in the UV). 



5 CONCLUSIONS 

Using GALEX UV and 2MASS J band photometry, we 
have investigated the relationship between UV-IR colours 
and spectroscopically derived stellar population parameters 
(age, metallicity and Q-abundance) for red sequence galaxies 
in local clusters. 

We select galaxies using strict emission criteria to avoid 
contamination from galaxies with very recent star formation. 
We analyse the NUV colour-magnitude relation (CMR) for 
our two samples of quiescent galaxies (920 in NFPS; 156 in 
SSC), and flnd rms dispersions of 0.37 and 0.30 mag (intrin- 
sic scatter of 0.36 and 0.29 mag) respectively. This is similar 
to previously reported values of ~0.5 mag and is an order 
of magnitude larger than the scatter in the optical CMR 
(~0.05 mag). 

We compared the NUV-J colour to the spectroscopic 
stellar population parameters for 87 galaxies in the SSC 
sample and found the following: 

• There is a significant NUV-J vs metallicity trend, with 
a slope of 1.27 ± 0.23 and an rms dispersion of 0.32 mag. 

• There is only a weak NUV-J vs age trend after the 
metallicity effect has been removed, and no correlation with 
a-abundance. 

• There is a large intrinsic scatter (~0.25 mag) in the 
NUV-J colour at fixed age and metallicity which cannot 
be easily accounted for with simple stellar populations. 

The unexpected blue colours of at least some objects, 
including an infiuential outlier, cannot be attributed to large 
radius contamination from other objects, and aperture bias 
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cannot account for the large scatter. Corrections for the UV 
upturn (UVX) phenomenon are relatively small (~0.2 mag) 
and are similar galaxy-to-galaxy, so do not reduce the in- 
trinsic scatter. 

We find that the large NUV-J intrinsic scatter could 
be attributed to galaxy 'frosting' by small (< 5 per cent) 
populations of either young stars or a low metallicity blue 
horizontal branch. 
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